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The lattice dynamics of LaFeAsOi-^Fz (a;=0, 0.1) and PrFeAsOi-y (y ~0.1) are investigated 
using inelastic x-ray scattering and ab-initio calculation. Measurements of powder samples 
provide an approximation to the phonon DOS, while dispersion is measured from a single 
crystal of PrFeAsOi-y. A model that agrees reasonably well with all of the data at room 
temperature is built from results of ab-initio calculations by reducing the strength of the Fe-As 
bond by 30%. 
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The recent discovery of high-temperature supercon- 
ductivity in iron-arsenic materials, 1-3 with T c exceeding 
50K in some compounds, 4,5 has led to an outpouring 
of work reminiscent of that following the discovery of 
the high T c of MgbV 6 In the latter case, the synergy of 
experiment and calculation quickly demonstrated that 
\ MgB2 was an exotic phonon-mediated superconductor. 
Despite the complexity of the new Fe-As compounds as 
compared to MgB2, calculations of phonon properties 
have appeared quickly, 7,8 and seem to support the as- 
sertion that these materials are not phonon mediated. 8,9 
However, these calculations require experimental verifi- 
cation. Here we show that measurements of the phonon 
density of states (DOS) in several samples, and phonon 
dispersion in a single crystal of PrFeAsOi_ y do not agree 
with LDA-based calculations, including those referenced 
above. However, the calculations can be brought into 
much closer agreement with the data by softening the Fe- 
As bond. This highlights the failure of the calculations 
to predict the properties of the most important bond in 
this system, and suggests that the results of those calcu- 
lations, especially regarding phonon behavior, should be 
carefully considered. 

It is early to survey the properties of the FeAs super- 
conductors, as information is very quickly evolving. How- 
ever, the available phonon calculations based on pseu- 
dopotential methods are essentially consistent with each 
other, 7, 8 and with our calculations here, and do not indi- 
cate strong electron-phonon coupling. 8 Calculations us- 
ing a rigid ion model also suggest small electron-phonon 
coupling. 9 The phonon DOS of a material provides a first, 
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global, look at phonon behavior. In this case, it allows us 
to develop a model of the dynamics based on the soften- 
ing of just the Fe-As bond, which agrees reasonably with 
the data. The agreement is confirmed by comparison with 
spectra measured from a single crystal of PrFeAsOi_ y . 

All measurements were made using ~1.5 meV res- 
olution at 21.75 keV at BL35XU 10 of SPring-8. The 
LaFeAsO and LaFeAsOn.gFn.i powder samples were syn- 
thesized as discussed in ref. 1 and the PrFeAsOi_ a was 
prepared using a high- pressure growth as in ref. 5. X-ray 
diffraction showed the polycrystalline samples to be pre- 
dominantly single phase, and a careful four-circle inves- 
tigation of the single crystal, 150x100x20 [im 3 , showed 
it to be a single grain, with a mosaic spread of about 1°. 
T c of the samples were found to be 27K for the powder 
of LaFeAsOo.gFo.i and 48K for PrFeAsOi_j,, consistent 
with refs. 1,5 and 36K for the single crystal. About 50 mg 
of each powder was gently pressed into a tube with a di- 
ameter of 3 mm to form a flat surface for reflection-based 
inelastic x-ray scattering measurements (they were con- 
firmed to be un-oriented after pressing) while the single 
crystal was mounted on a thin glass capillary and mea- 
sured in a Laue geometry. 

Calculations of the phonon spectra of LaFeAsO were 
carried out using the relaxed tetragonal P4/nmm struc- 
ture from VASP 11 calculated using a PAW method 12 
and the GGA. Phonons were then calculated via the 
PHONON package using a direct method, 13 in contrast 
to the perturbation methods of refs. 7, 8. In our calcu- 
lations, the energy cut-off for plane waves was 550 eV, 
the spacing of k-points was less than 0.1 A -1 and the 
convergence condition was that the total energy differ- 
ence be less than 1 /xeV. The relaxed structural parame- 
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Fig. 1. Calculated and measured densities of states. The usual 
phonon density of states from (a) Boeri et aZ., 8 (b) Singh et al.. 7 
and (c) this work, (d) The pseudo-DOS and (e) the effect of 
softening only the Fe-As bond (see discussion in text), (f)-(h) 
The pseudo-DOS extracted from measurements of three samples 
(see the Appendix for details on data processing). Vertical lines 
are at 13, 23, 30 and 35 meV, while the horizontal arrow shows 
the 5 meV shift between the room temperature measurements 
and the nominal calculation. The improved agreement with the 
soft Fe-As model is clear. 




ters were found to be a=4.0219 A and c=8.6215 A with 
z La =0. 14496 and ^=0.63817, similar to refs. 7,8. 

Figure 1 shows both the results of our calculations, 
and the data from several samples. The upper panels 
show the agreement of our calculated phonon DOS with 
previous work, and also the evolution of the DOS to 
the "pseudo-DOS" to be compared with experiment. By 
pseudo-DOS we mean the calculated one-phonon x-ray 
scattering integrated over the experimental range of mo- 
mentum transfers. In this case, we use a momentum full 
width of 3.7 nm _1 (corresponding to the acceptance of 
the analyzer array) centered at a momentum transfer 
of about 82 nm _1 - the precise value (fluctuation < 1 
run ) was sample dependent and was chosen to reduce 
the elastic contribution due to Bragg peaks. The DOS 
to pseudo-DOS transformation leaves the fine structure 
nearly unaffected, while changing the over-all weighting 
of the peaks. What is crucial for our discussion is that 
the triplet peak structure, with a broad peak at 13 meV 
and narrower ones at 23 and 35 meV remains robust, 
independent of which quantity is plotted. 

The lower half of Fig. 1 compares the extracted 
pseudo-DOS (see the appendix for details on data treat- 
ment) to the calculation. For all compounds, the first and 
second peaks agree reasonably with the calculation for 
LaFeAsO, however, the measured third peak is strongly 
softened, about 5 meV or 15%. Looking more closely, 
one also notes that the first peak is narrower experimen- 



Fig. 2. Selected mode polarizations and dispersion. Polarizations 
are shown for the points indicated on the dispersion plots, with 
modes at M and X shown with two unit cells to allow the full 
phasing to be visualized. Modes at V are generally either planar 
or c-axis polarized, while modes at M and X can be mixed. Modes 
at X, in particular, are complex and can involve circular motions. 
Red (blue) shading indicates Fe-As (oxygen) modes. The effect 
of the soft Fe-As bond is primarily to shift the Fe-As modes to 
lower energy, as compared to the original calculations. 

tally than calculated, while the second and third peaks 
are blurred out, and generally weaker than the calcu- 
lated ones. We also note the presence of a weak peak at 
about 39 meV in some of the measured spectra, and also 
seen, at slightly higher energy, in the calculation. Over- 
all, the clearest result is that the third peak in the DOS is 
softened. Considering the loss of definition of the peaks, 
however, one can speculate that the phonon modes be- 
tween 10 and 35 meV are somewhat blurred out and soft, 
either due to band broadening (increased dispersion), or, 
perhaps, lifetime effects. Our measured DOS results are 
consistent with ref. 14 (for E <20 meV) and unpublished 
work. 15, 16 

The dispersion and the contribution to the phonon 
modes near 35 meV calculated peak appears somewhat 
complex (see Fig. 2 and also discussion in ref. 8). The rel- 
evant band of modes consists primarily of in-plane FeAs 
motion at T, but the dispersion of these modes through 
the zone is complicated by (1) mixing with the c-axis 
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Fig. 3. Comparison of the measured dispersion for PrFeAsOi-^ 
(y ~0.1), T C =36K at room temperature (right panel) and calcu- 
lation (left). See text for discussion. 

polarized motion of the FeAs layers, and (2) the crossing 
of fast-dispersing oxygen modes. However, the peak in 
the pseudo-DOS is due nearly entirely to Fe-As motion. 
This can be demonstrated by softening only the Fe-As 
force constant matrix in the calculations. The resulting 
"soft Fe-As" pseudo-DOS calculation then agrees reason- 
ably well with the data: a reduction of the Fe-As nearest 
neighbor force constant matrix to 70% of the calculated 
value brings the calculated peak into line with the mea- 
sured one, and has only small effect on the rest of the 
spectrum. 

This model is in reasonable agreement with measure- 
ments of a small, 150x100x20 /im 3 , single crystal of 
PrFeAsOi_ y (y ~0.1), T C =36K. Figure 3 shows the dis- 
persion along r-X in the (300) Brillouin zone. The high 
energy region shows the softening of the LO modes, es- 
sentially as expected from our model, excepting near the 
transfer (£/=3.31, 3.42) of intensity from mode 3 in Fig. 
2 (an in-plane polarized LO mode) to mode 8 (a cir- 
cular mode in the a-c plane), where the data shows a 
more abrupt transfer than the calculation. The transi- 
tion point is sensitive to small changes in polarization, 
and could be affected by the oxygen nonstoichiometry. 
The slightly complex structures at lower energy come 
from exciting multiple nearby modes and are, partly, due 
to our finite momentum resolution [AQ= (0.067, 0.068, 
0.035) rlu]. However, interestingly, the broad structure 
at about 10 meV at (3.42, 0, 0.05) is actually in bet- 
ter agreement with the original calculations (not shown) 
and, will be the subject of future investigation. The re- 
maining small discrepancies in the lower-energy modes, 
given the changes in system (La^Pr), are not unreason- 
able. We note that available Raman data 17, 18 are in bet- 
ter agreement with the soft Fe-As model than the orig- 
inal calculation: the average energy difference between 
measurements and our calculation is reduced from 1.9 to 
0.9 meV (for three modes) 17 and from 2.0 to 0.9 meV 
(for five modes) 18 when the softening is included in the 
model. 

The softening of the phonons relative to the calcula- 
tions suggests that there is some character of the interac- 
tion that is not properly included in the calculations, and 



this interaction screens the phonons and reduces their 
energy. This is consistent with the tendency the LDA 
calculations to under-estimate the Fe-As bond distance 
(ie: 2.337 A in the present calculation as compared to 
2.407 A measured in ref. 19). In this context we note 
that another large discrepancy between the calculations 
and measured data regards the determination of the mag- 
netic moment on the iron atoms: the pseudopotential cal- 
culations 20 give consistently much larger iron moments 
(1 to 2 n B ) than measurements (~0. 2-0.4 ^ B ). 19 ' 21 ' 22 
This has been, improved by going to an all-electron cal- 
culation, 23 but bond-lengths were not reported. Also, a 
negative U 24 has been considered to account for this, but 
appears not to strongly affect the phonon energies. One 
pertinent fact, at least for the phonon question, regards 
the phase transition in the parent material LaFeAsO: it 
transforms from a tetragonal high temperature phase to 
an orthorhombic low temperature phase at about 145 
K, with magnetic order setting at ~T40K (see, e.g., refs. 
19 and 21). Carefully relaxed calculations suggest that 
the energy barrier between the two possible orthorhom- 
bic configurations at low-temperature is very small (~5 
meV), 25 well within phonon energy scales. Thus it is pos- 
sible that anharmonicity plays a significant role in deter- 
mining material properties. To investigate the effects of 
possible anharmonicity, we cooled the LaFeAsOo.gFo.i 
powder sample to 35K. We found the low-energy peak at 
13 meV hardened slightly, < 1 meV, the peak at 23 meV 
hardened about 1 meV, and the peak at 30 meV hard- 
ened about 1.5 meV, similar to refs. 15, 16. In order to 
get good agreement with calculation for the high-energy 
peak, we still have to soften the Fe-As force constant 
matrix by about 20%. While this does not rule out an 
anharmonic potential leading to softening of the Fe-As 
bond, it does clearly indicate that the strong softening 
we observe is predominantly intrinsic, not related to the 
room-temperature measurement. 

In summary, our data shows that one peak in the 
phonon DOS at room temperature is softened by about 
15% compared to LDA based calculations of LaFeAsO, 
and that this softening can be well modeled by reducing 
just the Fe-As nearest-neighbor force constant matrix by 
about 30% from the calculated value. This model is con- 
sistent with single crystal measurements. The softening 
of the peak in the pseudo-DOS is reduced to about 10% 
at 35K, consistent with about 20% reduction in the Fe- 
As force constant, but is still large. We also see a general 
blurring out of the spectral peaks that could be consis- 
tent either with increased dispersion of the modes (band- 
broadening) or possibly lifetime broadening, though the 
latter is not likely based on the single crystal work. 
These suggest rather strongly that the phonon system 
of this material is not yet understood and that perhaps 
one should not yet discard a possible phonon mediated 
mechanism on the basis of the presently available cal- 
culations. If one were to assume the only effect of the 
softening on the Eliashberg function, a 2 F, is to reduce 
some of the mode frequencies, then indeed it would make 
negligible change to the calculated T c , 8,15 however, the 
source of the softening is not clear, and, if it is due to 
coupling to the electronic system (so changing the mode 
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dependent electron-phonon coupling) then the effect on 
calculated T c could be larger. Notably, all samples inves- 
tigated, including the non-superconducting parent com- 
pound, show similar phonon DOS, suggesting that some 
other change in, e.g., the electronic structure, switches 
on superconductivity with doping. 
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Appendix: Data Analysis 

The extraction of the "pseudo-DOS" from the mea- 
sured scattering proceeds much has been done previously 
in neutron scattering. The procedure is exact in the case 
of pure incoherent scattering from a monoatomic har- 
monic solid at a well-defined momentum transfer, but 
approximate if there are significant other scattering pro- 
cesses, such as, in the present case, coherent scattering. 
A large range of integration over momentum space was 
chosen to attempt to maximize the effective incoherent 
contribution and an iterative fitting process was used to 
remove multiphonon contributions. This provides an ap- 
proximation to the generalized DOS (GDOS) which we 
call the pseudo-DOS in the paper. Early work along these 
lines can be found in ref. 26 while a recent application to 
inelastic x-ray scattering, essentially based on the anal- 
ysis for nuclear resonant scattering, is discussed in ref. 
27. This treatment preserves most of the main features 
of the spectra, while scaling them more nearly in agree- 
ment with the true one-phonon scattering (see Fig. A-l). 
One should note that the extracted spectra at low ener- 
gies are sensitive to the subtraction of the elastic peak, 
so that results below 8 meV should be taken as approx- 
imate. 
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Fig. A-l. Analysis progression for one data set at 300K. The top 
panel shows the measured data. The middle panel shows a scaled 
version, after subtraction of the elastic peak, that would approx- 
imate the GDOS in the limit of pure incoherent single-phonon 
scattering. The lowest panel shows the extracted pseudo-DOS 
after the multi-phonon contributions were removed by an itera- 
tive procedure. Note the solid line in the top panel is the back- 
calculated scattering from the pseudo-DOS, and its good agree- 
ment with the measurement demonstrates the self-consistency of 
the procedure. 
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